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Beschreibung 



Verfahren und Anordnung zur Ermittlung von Parametern eines 
technischen Systems 

Die Erfindung betrifft ein Verfahren sowie einer Anordnung zur 
Ermittlung von Parametern eines technischen Systems 

Bei einer Vielkanal-Obertragung und einem Vielkanal-Empf ang von 
Signalen tritt haufig eine Interferenz beispielsweise zwischen 
den Signalen/Bildern auf . Ein typisches Beispiel ist dabei eine 
Mischung eines Sprachsignals mit Rauschen, was bei der 
Telekommunikation und bei Videokonf erenzen ein grofJes Problem 
darstellen kann. Die vorliegende Erfindung betrifft daher das 
Gebiet der Signaltrennung, um beispielsweise ein ursprungliches 
SprachsiLgnal wiederzugewinnen. 

Typische bekannte Techniken zur Trennung der Quellensignale 
ausgehend von Mischsignalen basieren auf einer zeitlichen 
Mittelung oder einer Filterung der Signale. Dies hat indessen 
Nachteile beziiglich des Rechenauf wands zur' Folge. 

Es sind auch Verfahren auf der Grundlage der sogenannten Blind 
Channel Equalization (Signalent zerrung ohne Vorkenntnisse iiber 
den Ubertragungskanal ) bekannt, aber diese Verfahren benotigen 
immer eine gewisse Kenntnis iiber die Quellsignale, wie 
beispielsweise eine Kenntnis iiber ihre statistische Verteilung. 

Das Problem der Signaltrennung tritt beispielsweise auch auf, 
wenn zwei Sprecher in zwei entfernt von ihnen stehende 
Mikrofone sprechen, so daB jedes Mikrofon ein Gemisch der von 
den zwei Sprechern gesprochenen Signale aufnimmt. Somit besteht 
das Problem, das Signalgemisch, also eine Menge uberlagerter 

-Eingangss.ignale— wi-eder— zu— trennen- Au^-^ L— Moiged-e-y-,— H-.-G 

Schuster, ^Separation of a Mixture of Independent Signals using 
Time-Delayed Correlations", Phys. Ref. Lett. 72, 3634 (1994) 
ist dabei das folgende Verfahren bekannt: Das Problem, n 
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iiberlagerte und korrelierte Quellensignale (Eingangssignale) zu 
trennen und gleichzeitige die Mischungskoef f izienten der 
Quellenstarken zu bestimmen, lafit sich auf ein Eigenwertproblem 
reduzieren, bei dem simultan zwei symmetrische n x n Matrizen 
diagonalisiert werden miissen. Die Matrixelemente sind meftbare 
zeitverzogerte Korrelationsf unktionen. Die Losung- des 
Eigenwertproblems kann durch ein neuronales Netz erfolgen, 
wobei die Lernregel fur die lateralen inhibitorischen 
Wechselwirkungen zwischen den Neuronen durch eine- Liapunov- 
Funktion bestimmt werden kann, deren Minima die (entarteten) 
Losungen des Problems liefern. 

Dieses Verfahren wurde auch bereits auf akustische 
Eingangssignale angewandt (siehe F. Ehlers, H. G. Schuster, 
„Blind Separation of convolutive mixtures and an application in 
automatic speech recognition" IEEE Trans. Signal Proc. (1997) . 

Aus der DE 195 31 388 CI ist ein Signaltrennungsverf ahren und 
eine Signaltrennungseinrichtung fur nichtlineare Mischungen 
unbekannter Signale (Blind Channel) bekannt, das schematisch in 
Fig. 3 dargestellt ist. 

Dieses deutsche Patent behandelt die Separierung eines 
Signalgemisches, bestehend aus der nichtlinearen Uberlagerung 
von M unbekannten Quellsignalen XI, X2, wobei N (N > M) 
unterschiedliche Mischungen der M Quellsignale XI, X2 inklusive 1 
eines eventuellen Storsignals einer Signalauswerteeinrichtung 
zugefuhrt werden, die durch eine statistische Analyse der 
Signale die nichtlinearen Obertragungsf aktoren bestimmt und mit 
diesen errechneten Faktoren die Mischung riickgangig macht, so 
daft die N Ausgange der Signaltrennungseinrichtung moglichst 
naherungsweise die M Quellsignale ohne Oberlagerungen 
enthalten. Dadurch wird eine Behandlung nichtlinearer Gemische 
moglich, wobei nichtlinear bedeutet, dafi die Quellsignale XI, 

"X 2" ~du r c h~~ei n u nbe k a n n t e s n i^hirl~iTr^axes~" S y s ; t"e"m"~G~g"e"m i~s"c h t — 

werden. Das unbekannte System G wird durch eine sogenannte 
Volterra-Reihe beschrieben, und die Signaltrennungseinrichtung 
G-l bestimmt die Koef f izienten der Volterra-Reihe. Mit deren 
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Kenntnis ist eine Entmischung des Signalgemischs moglich. 
Auflerdem konnen die Koef f izienten zu weiteren Analysen zur 
Orts- oder Geschwindigkeitsbestimmung der Signalquellen benutzt 
werden. 

Das aus dieser Druckschrift bekannte Verfahren besteht dabei im 
wesentlichen aus zwei Schritten: 

• Erstens werden die durch den wahlbaren Grad der* 
Nichtlinearitat bei der Mischung eindeutig bestimmten 
nichtlinearen Gleichungen fur ein gleitendes Zeitf enster 
gelost und die Losungen werden liber die Zeit gemittelt. Diese 
zeitliche Mittelung stellt einen Hauptnachteil dieser 
bekannten Technik dar, da sie den Rechenauf wand und 
gleichzeitig die Zeitdauer fur die Berechnung erhoht. 

• Zweitens wird aus einer genugend grolien Anzahl von 
unterschiedlichen Kumulanten der geschatzten Ausgangssignale 
gebildetes Potential minimiert , wobei die zur Berechnung des 
Potentials notigen Werte aus einem gleitenden Zeitfenster 
einer wahlbaren Lange stammen. Dabei ist vorausgesetzt, daft 
sich das Mischungssystem so langsam andert, daft diese 
Anderung bei der Berechnung der gesuchten Mischungsf aktoren 
vernachlassigt werden kann. Gemafi diesem deutschen Patent 
wird zur Ausfuhrung des zweiten genannten Schrittes eine 
Kostenf unktion konstruiert, die minimiert wird. Wenn das 
globale Minimum erreicht ist, hat man die optimalen Werte, in 
diesem Fall der Ubertragungsf aktoren, gefunden. 

Hinsichtlich des zeitlichen Aufwands sowie des Rechenauf wands 
ist das in der DE 195 31 388 CI beschriebene Verfahren aufgrund 
der Ausfuhrung einer zeitlichen Mittelung am Abschluft des 
ersten, oben genannten Verf ahrensschrittes nachteilig. 

~Dre- vorl-regende—Exfxndung - hair - daher— zur~ Aufgabe- — e±n~ Vexfa+rren 
und eine Anordnung bereit zustellen, die die Trennung 
liberlagerter , statistisch voneinander unabhangiger akustischer 
Signale mit verringertem Rechenauf wand ermoglicht. 



GR 98 P 2958 



Diese Aufgabe wire! durch das Verfahren sowie durch die 
Anordnung mit den Merkmalen gemaft den unabhangigen Anspruchen 
gelost. 

Bei einem Verfahren zur Ermittlung von Parametern. eines 
technischen Systems, mit dem Ausgangssignale aus einer Menge 
iiberlagerter , statistisch voneinander unabhangiger 
Eingangssignale ermittelt werden konnen, werden die Parameter 
derart ermittelt, daft die statistische Unabhangigkeit der 
Ausgangssignale maximiert wird. 

Eine Anordnung zur Ermittlung von Parametern eines technischen 
Systems, mit dem Ausgangssignale aus einer Menge iiberlagerter , 
statistisch voneinander unabhangiger Eingangssignale ermittelt 
werden konnen, weist einen Prozessor auf, der derart 
eingerichtet ist, daft die Parameter derart ermittelt werden 
konnen, daft die statistische Unabhangigkeit der Ausgangssignale 
maximiert wird - 

Vorteilhafte Weiterbildungen der Erf indung ergeben sich aus den 
abhangigen Anspruchen. 

Die Parameter werden bevorzugt in einem iterativen Verfahren 
ermittelt . 

In einer weiteren Ausgestaltung sind die Parameter Elemente 
einer Entmischmatrix, mit der die Menge der iiberlagerten 
Eingangssignale multipliziert oder auch gefaltet wird, wodurch 
die Ausgangssignale gebildet werden. 

Die Optimierung der Parameter der Entmischmatrix wird bevorzugt 
durch die folgenden Schritte er halten : 

- Wiederholung einer zeitverzogerten Dekorrelationsberechnung 
zur Ermittlung der Eigenwerte der Entmischmatrix, 
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- Ermittlung der Eigenwerte der Entmischmatrix, fur die 
Kreuzkorrelationen einen minimalen Wert annehmen, und 

- Ausfiihrung einer Kumulantenminimierung, wobei als Startwerte 
fur die Kumulantenminimierung die im vorherigen Schritt 
ermittelten Eigenwerte verwendet werden. 

Die Kumulantenminimierung kann beispielsweise durch Trainieren 
eines neuronalen Netzes, aber auch jegliche andere bekannte 
Minimierungstechnik, wie beispielsweise Gradientenabstieg oder 
Monte-Carlo-Simulationen verwendet werden. 

In einer Weiterbildung wird bei der Optimierung der Parameter 
der Entmischmatrix wenigstens ein Diagonalparameter der 
Entmischmatrix auf einen vorgegebenen Wert gesetzt, womit , um 
eine Stabilitat des Minimierungsvorgangs hin zu einem globalen 
Minimum zu gewahrleisten . 

Die Entmischmatrix wird bevorzugt auf eine finite Impulsantwort 
begrenzt, d.h. es wird ein FIR-Filter (Finite ]jmpulse Response) 
eingesetzt zur Bildung der einzelnen Komponenten der 
Entmischmatrix. Der FIR-Filter kann sowohlein kausales FIR- 
Filter oder auch ein nicht-kausales FIR-Filter sein. 

Ferner wird die Entmischmatrix wahrend der 

Kumulantenminimierung bevorzugt durch Projektion in einen 
Einheitskreis stabilisiert . 

Die Weiterbildungen gelten sowohl fur das Verfahren als auch 
fur die Anordnung, wobei bei der Anordnung jeweils der 
Prozessor derart eingerichtet ist, daft der entsprechende 
Verf ahrensschritt durchfuhrbar ist oder durchgefuhrt wird. 

D i e_E r f i ndun g-_s.O-w±e__dex_e n__We i.t.e r bi 1 dun g.e n_konn e n__v.orXjei.lha.f_t 

eingesetzt werden zur Trennung uberlagerter , statistisch 
voneinander unabhangiger Eingangssignale, insbesondere 
akustischer Eingangssignale . 
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Das Verfahren sowie die Anordnung sind fur eine beliebige 
Anzahl von Eingangssignalen anwendbar. 

Weitere Vorteile, Merkmale und Eigenschaf ten der vorliegenden 
Erfindung werden nunmehr bezugnehmend auf die beiliegenden 
Figuren der Zeichnungen anhand eines Ausf uhrungsbeispiels" nahei 
erlautert. 

Fig. 1 zeigt die Anwendung eines Systems zur Trennung 
uberlagerter, statistisch voneinander unabhangiger 
akustischer Signale gemaii dem Ausf iihrungsbeispiel, 

Fig. 2 zeigt eine symbolische Darstellung des Systems von 
Fig. 1, und 

Fig. 3 zeigt eine aus dem Stand der Technik (DE 195 31 388 
CI) bekannte Signaltrennungseinrichtung fur nichtlineare 
Mischungen unbekannter Signale. 

Zur Wiedergewinnung des ursprunglichen Sprachsignals wird aus 
einer Mischung von Signalen die statistische Unabhangigkeit 
zwischen den Quellensignalen (ursprungliches Sprachsignal und 
Rauschen) , im weiteren auch als Eingangssignale bezeichnet, 
ausgenutzt und der inverse Vorgang des dynamischen Systems, de 
die Mischung der Signale ergeben hat, wird im wesentlichen 
naherungsweise trainiert (gelernt) . Zwei verschiedene 
Mischungen des Sprachsignals bzw. des Rauschsignals werden 
beispielsweise durch zwei Mikrofone 1, 2 (vgl. Fig.l) erhalten, 
die voneinander beabstandet sind und/oder in entgegengeset zten 
Richtungen ausgerichtet sind. Bei dem Verfahren wird die 
sogenannte zeitverzogerte Dekorrelationstechnik (TDD, time 
delayed decorrelation) verwendet, urn die Lernphase zu 
initiieren, d.h. urn Startwerte fur die Lernphase zu ermitteln 



und vorzugeben, wodurch der Berechnungsauf wand fur eine im ~~ 
weiteren beschriebenen Kumulantenminimierung verringert -werden 
kann und die Gefahr lokaler Minima verringert werden kann. 
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Fig.l zeigt zwei Mikrofone 1, 2, die ein erstes Eingangssignal 
Zl(t) und ein zweites Eingangssignal Z2 (t) aufnehmen. Diese 
Eingangssignale Zl(t) und Z2(t) konnen untereinander wiederum 
jeweils mit Rauschen vermischt sein, was durch eine Mischmatrix 
S (siehe Bezugszeichen 3) symbolisch in Fig.l dargestellt ist. 
Nach dem Empfang bzw. der Obertragung wird eine Menge XI (t) und 
X2(t) iiberlagerter , statistisch voneinander unabhangiger 
Eingangssignale Zl(t) und Z2(t) erhalten. Diese Signale werden 
in eine Berechnungseinheit 4 eingegeben, in der im wesentlichen 
zwei Schritte ausgeftihrt werden, die symbolisch durch eine 
Berechnungseinheit B (Bezugszeichen 6) fur den ersten Schritt 
sowie ein neuronales Netzwerk 5 fur den zweiten Schritt 
dargestellt ist. 

Durch die Berechnungseinheit 4 werden zwei Ausgangssignale 
Yl(t) bzw. Y2(t) ermittelt, die bei optimaler Einstelllung der 
Parameter in der Berechnungseinheit 4 naherungsweise gleich den 
Eingangssignalen Zl(t) bzw. Z2(t) sind. Mit anderen Worten, bei 
optimaler Einstellung der Parameter der verwendeten Matrizen in 
der Berechnungseinheit 4 erfolgt durch diese Berechnungseinheit 
4 im wesentlichen der inverse Vorgang zu dem dynamischen 
Mischvorgang, der durch die Matrix S (Bezugszeichen 3) 
symbolisch dargestellt ist. Das Ausf iihrungsbeispiel beschaftigt 
sich mit dem Optimierungsvorgang der Einstellung der Parameter 
der Entmischmatrix . 

Es werden die Parameter der Matrizen in der Berechnungseinheit 
4 derart optimiert, daft die statistische Unabhangigkeit 
zwischen den durch den Matrizenvorgang in der 
Berechnungseinheit 4 gewonnenen Ausgangssignale Yl(t), Y2(t) 
maximiert wird. Zu diesem Zweck werden die Ausgangssignale 
Yl(t) bzw. Y2 (t) zu der Berechnungseinheit 4 zuruckgef uhrt (s. 
Ruckf uhrschleif en 7 bzw. 8) . Durch ein iteratives Verfahren 
wird ermittelt, ob sich die statistische Unabhangigkeit der 
Ausgangssrgnal'e - Yl~(t~) bzw~ Y2~(t ) im~Ve~rgre~rchn~zu~ dem~ vorherigen 
Schritt der iteration erhoht hat (und somit die Iteration die 
„richtige" Richtung in Richtung des globalen Minimums einer im 
weiteren beschriebenen Kostenf unktion einnimmt) oder nicht. 
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Fig. 2 zeigt eine mathematische Darstellung des Schemas von 
Fig.l, wobei der Mischvorgang 3 durch eine Matrix S (q) 
mathematisch beschrieben werden kann und der Entmischvorgang, 
der durch die Berechnungseinheit 4 erfolgen soil, durch eine 
Entmischmatrix M(q) symbolisiert wird. 

In Fig. 2 ist somit das Problem der Trennung einer sogenannten 
Multi-Channel Blind Source (Vielf ach-Kanal-Quelle - ohne a- 
priori-Kenntnis) in zwei Dimensionen dargestellt. Dabei ist 
angenommen, dali das Mischsystem S (q) , wobei q fur eine 
Einheitsverzogerung ste'ht, stabil ist und gleichzeitig auch 
eine stabile Invertierung aufweist, d.h. dali es ein 
Minimalphasensystem ist. Daruber hinaus ist angenommen, dali di 
Eingangssignale Zl(t) und Z2(t) (beispielsweise ein Sprach- 
bzw. ein Rauschsignal) statistisch voneinander unabhangig sind 
und nicht gaussformig verteilt sind. Die Menge XI (t) und X2(t) 
der uberlagerten Eingangssignale Zl(t) und Z2 (t) sind 
Eingangssignale in ein Entmischsystem mit einer Entmischmatrix 
M(q), deren Parameter (Matrixelemente) auf eine Maximierung de 
statistischen Unabhangigkeit zwischen den Ausgangssignalen 
Yl(t) und Y2(t) trainiert werden. Unter „Trainieren" ist dabei 
der gut bekannte Lernvorgang beispielsweise eines neuronalen 
Netzes bekannt, das als ein Beispiel fur eine Technik genannt 
sein soil, die statistische Unabhangigkeit zu maximieren. Dies 
erfolgt durch Minimierung einer im weiteren beschriebenen 
Kostenf unktion J(M) . 

Es wird eine Kumulanten-Kostenf unktion gebildet, die die 
Diagonalkumulantenelementen der Kumulantenordnung 2-4 
minimiert : 



Dcum « J(M) - 



4 

i = l nondiag 



C^nondiag 



Folgende Gesichtspunkte der dynamischen Mischung durch die 
Mischmatrix S(q) sind dabei zu berucksichtigen : 
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Stabilitat des Entmischsystems : 
Dies wird erreicht, wenn M(q) auf eine finite 
Impulsantantwort (FIR-Filter) beschrankt wird. Die 
Stabilitat des FIR-Systems M(q) kann daruber hinaus auch 
dadurch erhalten werden, daB wahrend der Lernphase eine 
Projektion in den Einheitskreis erfolgt Eine 
moglicherweise vorliegende Nichtkausalitat der 
Invertierung von S(q) kann durch eine geeignete 
Zeitverschiebung ( Verzogerung) des Eingangssignals X(t) 
kompensiert werden. 

Eindeutigkeit der getrennten Signale Y(t): 
Fur den Fall statischer Mischungen werden die 
ursprunglichen Quellsignale durch eine Skalierung 
wiedergewonnen. Fur den Fall einer dynamischen Entmischung 
ist die Gefahr einer Nichteindeutigkeit der getrennten 
Signale Y(t) sogar noch groBer. Es ist of f ensichtlich, daft 
fur den Fall, daB Yl(t) und Y2(t) statistisch voneinander 
unabhangig sind, auch jegliche linear gefilterte 
Modifikation dieser Signale immer noch statistisch 
unabhangig sind. Daher ist eine zusatzliche Information 
notwendig, urn die inharente Nichteindeutigkeit des 
Problems zu verringern. 

Gaussverf ormung der Daten: 

Algorithmen auf Kumulantenbasis fur eine statische Blind 
Source-Trennung minimieren bzw. eliminieren in effektiver 
Weise Diagonalkumulanten hoherer Ordnung entsprechend den 
Ausgangssignalen Y(t). Andererseits fiihrt eine lineare 
Filterung zu einer Gaussverteilungs-Verf ormung der Daten, 
bei denen die Kumulanten hoherer Ordnung in Richtung 0 
gehen. Dies kann daher zu Randlosungen fiihren, in denen 
die Kostenf unktion ein lokales Minimum erreicht, wobei die 
erwiinschte - ei"gentiL~i"ch"e~Trennung (~gl~obaTes "Minimum) ni~cht 
erfolgt. Urn diesen ungewunschten Fall zu vermeiden, wird 
die Struktur der Entmisch-Transf erf unktion 

(Entmischmatrix) M(q) einigen Beschrankungen unterworfen. 
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Urn die obengenannten Probleme zu umgehen, wird ein Ansatz 
gewahlt, dali wenigstens eines (oder auch alle) der 
Diagonalelemente auf den Einheitswert gesetzt werden: 



Mll(q) = 1 



und 

M22(q) = 1. 

Diese Annahme ist exakt, wenn die Mischelemente Sll (q) und/ode 
S22(q) ebenfalls einen Einheitswert „1" aufweisen. Sonst sei 
angenommen, daft Sll (q) und/oder S22 (q) stabile Invertierungen 
aufweisen, was die Skalierung der Diagonalelemente von M(q) au 
den Einheitswert erlaubt. Dieser Ansatz verringert wesentlich 
die Nichteindeutigkeit einer Losung und vermeidet dariiber 
hinaus in effektiver Weise die Gefahr einer ubermafiigen 
Gaussverteilungs-Verf ormung der Ausgangssignale . Auch wenn auf 
den ersten Blick die Beschrankung der Diagonalelemente von 
M(q), wie oben ausgefuhrt, sehr restriktiv wirkt, ist diese 
Annahme in der praktischen Anwendung in der Regel erfullt. Ein 
typisches Beispiel ist die Entrauschung von Sprachsignalen auf 
der Grundlage einer Aufzeichnung mit zwei Mikrofonen, wobei die 
Mikrofone voneinander raumlich getrennt sind oder ein Mikrofon 
in Richtung des Sprechers gerichtet ist, wahrend das andere 
Mikrofon in der umgekehrten Richtung gerichtet ist, so dali das 
zweite, von dem Sprecher abgewandte Signal im wesentlichen nur 
ein Rauschsignal auf nimmt . 

Der Kumulantenansatz basiert auf einer direkten Ermittlung der 
Diagonalkumulanten, wie sie in dem eingangs genannten Artikel 
von F. Ehlers und H. Schuster, „Blind Separation of convolutive 
mixtures and an application in automatic speech recognition" 

IEEE"Trans - STgnaT~Proc; (T9"9T) , ausgefQHrf - i"sf was aber 

inharent den Nachteil aufweist, dali die numerische Losung sehr 
aufwendig ist. Daher wird eine geeignete Initialisierung dieses 
Minimierungsverf ahrens verwendet. Urn Startwerte fur das 
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Minimierungsverf ahren zu ermitteln, wird gemaft der vorliegenden 
Erfindung die Technik der zeitverzogerten Dekorrelation (TDD) 
zur gleichzeitigen Dekorrelierung zweier verschiedener 
Zeitverzogerungen angewendet, wobei diese TDD-Technik sich auf 
ein geeignetes Matrix-Eigenwertproblem zuriickfiihren laftt. Wie 
bereits gesagt, wird diese TDD-Technik gemafi der vorliegenden 
Erfindung zur Initiierung des Diagonal (Kreuzkorrelations ) - 
Kumulantenminimierungsproblems verwendet . 

Zusammenf assend lafit sich das Verfahren in die folgenden beiden 
Schritte unterteilen: 

1. Wiederholung des TDD-Verf ahrens auf Grundlage des 
Eigenwertproblems im Frequenzbereich fur verschiedene 
Verzogerungspaare und Ermittlung der Losung, bei der die 
Kreuzkorrelationsterme einen minimalen Wert aufweisen. 

2. Initiierung (Start) der Diagonal-Kumulantenminimierung auf 
Grundlage der im oben genannten Schritt ermittelten 
Startwerte ( FIR-Parameter ) . 

Im folgenden sollen noch einige Haupteigenschaf ten und Vorteile 
noch einmal zusammengef aftt werden: 




Es ist kein a-priori-Wissen der Signaleigenschaf ten 
notwendig, mit der Ausnahme, daft eine statistische 
Unabhangigkeit gefordert wird. 

Die Stabilitat des dynamischen Entmischsystems wird durch 
die Modulierung seiner Komponenten als FIR-Filter 
gewahrleistet • 

Eine ubermaftige Gaussverteilung-Verf ormung wird durch den 
Ansatz vermieden, daft wenigstens eines der Elemente der 
~M*:Ts7?ht'r"a"n^^ (Entrnxschmatrrx) auf~~den 

Einheitswert gesetzt wird bzw. auf den Einheitswert 
skaliert werden kann, und 
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da der Kumulanten-Minimierungsschritt (Schritt 2) einen 
hohen Rechenauf wand erfordert, wird der Lernalgorithmus 
beispielsweise eines neuronalen Netzes durch das TDD- 
Verfahren initialisiert . 
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Im weiteren ist ein Programm in Matlab, Version 4 oder Version 
angegeben, mit dem das oben beschriebene Ausf uhrungsbeispiel 
auf einem Rechner realisiert werden kann: 



function [cost, outl, out2] =cumulant_ costFIRa2 (par, input, pi, pll, 
p2,p22,a3,a4) ; 

% [cost, outl, out2] =cumulant_costFIRa2 (par, input , pi, pll, p2, p22, 
a3, a4 ) ; 

% cumulant cost 

% FIR representation used 

% filter function used in both directions (non_causal) 

[np, mp] =size (par) ; 
f irl=par (1 :pl) ; 
firll=par ( 1+pl : pl+pll ) ; 
fir2=par ( 1+pl+pll : pl+pll+p2 ) ; 
f ir22=par (pl+pll+p2+l :mp) ; 

den=l; %FIR only 

out 1= [input ( : ; 1) -filter (fir 1, den, input ( : , 2) ) -f lipud ( filter 
([0 firll ], [den] , f lipud (input (: ,2) )))] ; 

%/std(input (1, : ) ) ; %dlsim 

%f ilter 

out2= [input ( : , 2) -f ilter ( fir2 , den, input ( : , 1) ) -f lipud ( filter 
([0 fir22 ] , [den] , f lipud (input ( : , 1) ) ) ) ] ;%/std (input ( : , 2) ) ; 

%dlsim %f ilter 

out=[outl out2]; 

%outl=outl/std (outl ) ; % this scaling was not needed in examples 
in SIP98 paper 
%out2=out2/std (out2) ; 

Ld=0 ; — % ~nuiTib"e"r~ of"~^ de~rays~rn—ca~l-cura ting" the cr os s-cor rel"at±on~ 

cost3=0; 

cost4=0; 
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costALLl=[] ; 
costALL2-[] ; 
ol2=outl. *out2; 
cost2=mean (ol2) A 2; 

oll2=outl. *outl. *out2; 
ol22=outl. *out2. *out2; 

if a 3 ==1 

cost3=[mean(ol22) ] A 2 + [mean (ol22 ) ] A 2; 

end 

if a4 ==1 

cost 4= [mean (oll2 . *outl) -3*mean (outl . A 2) *mean (ol2) ] 
[mean( (outl. A 2) .*(out2. A 2) ) -2*mean (ol2 ) A 2- 
mean (outl. A 2) *mean (out2 . A 2)] A 2+. . . 
[mean (o!22 . *out2) -3*mean (out2 . A 2) *mean (o!2) ] A 2; 



end 



%cum4a= [cum4x (outl, outl, outl, outl) ] A 2; 
%cum4b= [cum4x (out2 , out2 , out2 , out2 ) ] A 2; 



cost=cost2+a3*cost3+a4*cost4 ; 



%-cum4a-cum4b; 
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Patentanspruche 

1. Verfahren zur Ermittlung von Parametern eines technischen 
Systems, mit dem Ausgangssignale aus einer Menge iiberlagerter, 
statistisch voneinander unabhangiger Eingangssignale ermittelt 
werden konnen, bei dem die Parameter derart ermittelt werden, 
daft die statistische Unabhangigkeit der Ausgangssignale 
maximiert wird. 

2. Verfahren nach Anspruch 1, 

bei dem die Parameter in einem iterativen Verfahren ermittelt 
werden. 

3. Verfahren nach Anspruch 1 oder 2, 

bei dem die Parameter Elemente einer Entmischmatrix sind, mit 
der die Menge der uberlagerten Eingangssignale multipliziert 
wird, wodurch die Ausgangssignale gebildet werden. 

4. Verfahren nach Anspruch 3, 

bei dem die Optimierung der Parameter der Entmischmatrix durch 
die folgenden Schritte erhalten wird: 

- Wiederholung einer zeitverzogerten Dekorrelationsberechnung 
(6) zur Ermittlung der Eigenwerte der Entmischmatrix, 

- Ermittlung der Eigenwerte der Entmischmatrix, fur die 
Kreuzkorrelationen einen minimalen Wert annehmen, und 

- Ausfuhrung einer Kumulantenminimierung (5) , wobei als 
Startwerte fur die Kumulantenminimierung die im vorherigen 
Schritt ermittelten Eigenwerte verwendet werden. 

5. Verfahren nach Anspruch 4, 

bei dem die Kumulantenminimierung durch Trainieren eines 
neuronalen Netzes (5) erfolgt. 



6. Verfahren nach einem der Anspruche 3 bis 5, 
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bei dem bei der Optimierung der Parameter der Entmischmatrix 
wenigstens ein Diagonalparameter der Entmischmatrix auf einen 
vorgegebenen Wert gesetzt wird. 

7. Verfahren nach einem der Anspruche 3 bis 6, 

bei dem die Entmischmatrix auf eine finite Impulsantwort 
begrenzt wird. 

8. Verfahren nach einem der Anspruche 3 bis 7, 

bei dem die Entmischmatrix wahrend der Kumulantenminimierung 
(5) durch Projektion in einen Einheitskreis stabilisiert wird. 

9. Verfahren nach einem der Anspruche 1 bis 8, eingesetzt zur 
Trennung uberlagerter , statistisch voneinander unabhangiger 
Eingangssignale . 

10. Verfahren nach einem der Anspruche 1 bis 8, eingesetzt zur 
Trennung uberlagerter, statistisch voneinander unabhangiger 
akustischer Eingangssignale . 

11. Anordnung zur Ermittlung von Parametern eines technischen 
Systems, mit dem Ausgangssignale aus einer Menge uberlagerter, 
statistisch voneinander unabhangiger Eingangssignale ermittelt 
werden konnen, mit einem Prozessor, der derart eingerichtet 
ist, daft die Parameter derart ermittelt werden konnen, dafl die 
statistische Unabhangigkeit der Ausgangssignale maximiert wird. 

12. Anordnung nach Anspruch 11, 

bei der der Prozessor derart eingerichtet ist, dafi die 
Parameter in einem iterativen Verfahren ermittelt werden. 

13. Anordnung nach Anspruch 11 oder 12, 

"bei~~der der Prozessor derart eingerichtet ist, daft die 
Parameter Elemente einer Entmischmatrix sind, mit der die Menge 
der uberlagerten Eingangssignale multipliziert wird, wodurch 
die Ausgangssignale gebildet werden. 
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14. Anordnung nach Anspruch 13, 

bei der der Prozessor derart eingerichtet ist, daft die 
Optimierung der Parameter der Entmischmatrix durch die 
folgenden Schritte erhalten wird: 

- Wiederholung einer zeitverzogerten Dekorrelationsberechnung 
(6) zur Ermittlung der Eigenwerte der Entmischmatrix, 

- Ermittlung der Eigenwerte der Entmischmatrix, fur die 
Kreuzkorrelationen einen minimalen Wert annehmen, und 

- Ausfiihrung einer Kumulantenminimierung (5), wobei als 
Startwerte fur die Kumulantenminimierung die im vorherigen 
Schritt ermittelten Eigenwerte verwendet werden. 

15. Anordnung nach Anspruch 14, 

bei der der Prozessor derart eingerichtet ist, daft die 
Kumulantenminimierung durch Trainieren eines neuronalen Netzes 
(5) erfolgt . 

16. Anordnung nach einem der Anspriiche 13 bis 15, 

bei der der Prozessor derart eingerichtet ist, daft bei der 
Optimierung der Parameter der Entmischmatrix wenigstens ein 
Diagonalparameter der Entmischmatrix auf einen vorgegebenen 
Wert gesetzt wird. 

17. Anordnung nach einem der Anspriiche 13 bis 16, 

bei der der Prozessor derart eingerichtet ist, daft die 
Entmischmatrix auf eine finite Impulsantwort begrenzt wird. 

18. Anordnung nach einem der Anspriiche 13 bis 17, 

bei der der Prozessor derart eingerichtet ist, daft die 
Entmischmatrix wahrend der Kumulantenminimierung (5) durch 
Projektion in einen Einheitskreis stabilisiert wird. 



19. Anordnung nach einem der Anspriiche 11 bis 18, eingesetzt 
zur Trennung uberlagerter , statistisch voneinander unabhangiger 
Eingangssignale . 
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20. Anordnung nach einem der Anspruche 11 bis 18, eingesetzt 
zur Trennung iiberlagerter , statistisch voneinander unabhangige 
akustischer Eingangssignale . 
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Zusammenf assung 

Verfahren und Anordnung zur Ermittlung von Par ame tern eines 
technischen Sy s terns 

Es werden Parameter eines technischen Systems, mit dem 
Ausgangssignale aus einer Menge uberlagerter , statistisch 
voneinander unabhangiger Eingangssignale ermittelt werden 
konnen, bestimmt. Die Parameter werden derart ermittelt, dafi 
die statistische Unabhangigkeit der Ausgangssignale maximiert 
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